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Abstract. A eoneept of “eritioal” simplifieation was proposed by Yablonsky and Lazman in 1996 
IfTSlI for the oxidation of earbon monoxide over a platinum eatalyst using a Langmuir-Hinshelwood 
meehanism. The main observation was a simplifieation of the meehanism at ignition and extinetion 
points. The eritieal simplifieation is an example of a mueh more general phenomenon that we eall 
a bifurcational parametric simplification. Ignition and extinetion points are points of equilibrium 
multiplieity bifureations, i.e., they are points of a eorresponding bifureation set for parameters. 

Any bifureation produees a dependenee between system parameters. This is a mathematieal 
explanation and/or justifieation of the “parametrie simplifieation”. It leads us to a eonjeeture that 
“maximal bifureational parametrie simplifieation” eorresponds to the “maximal bifureation eom- 
plexity.” 

This eonjeeture ean have praetieal applieations for experimental study, because at points of 
“maximal bifurcation complexity” the number of independent system parameters is minimal and 
all other parameters can be evaluated analytically or numerically. 

We illustrate this method by the case of the simplest possible bifurcation, that is a multiplicity 
bifurcation of equilibrium and we apply this analysis to the Langmuir mechanism. Our analytical 
study is based on a coordinate-free version of the method of invariant manifolds (proposed recently 
0). As a result we obtain a more accurate description of the “critical (parametric) simplifications.” 

With the help of the “bifurcational parametric simplification” kinetic mechanisms and reaction 
rate parameters may be readily identified from a monoparametric experiment (reaction rate vs. 
reaction parameter). 
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1. Introduction and methodology 

By a standard linguistic definition, a bifurcation is a place where something divides into two 
branches. Typical examples are: bifurcation of aorta is the region in which the abdominal aorta 
bifurcates into the left and right branches; river bifurcation, the forking of a river into its distribu¬ 
taries; bifurcation of trachea, etc... For all these examples, bifurcations are places where the main 
stream is divided into two or more streams. This dynamical interpretation corresponds to geometric 
changes of objects. 

This simplistic meaning of bifurcation can be extended. Roughly speaking, a bifurcation is a 
region of qualitative changes in the system’s dynamical behavior. It can be induced by changes in 
geometry, but also by other causes. 

The notion of “bifurcation” in mathematics was introduced by Henri Poincare lfT3ll at the end 
of nineteen century for changes in the multuplicity of equalibrium systems of ordinary differential 
equations (ODE). Nowadays, this notion is used for systems of partial differential equations (PDF) 
and all other types of continuous and discrete mathematical models. Any mathematical model, 
except for a set of variables that describes its dynamical behavior, also has a set of parameters that 
are “permanent” comparatively with variables. In reality these parameters are control parameters 
and the system’s asymptotic dynamics (for example, the equilibrium multiplicity) depends on its 
values. 

Roughly speaking, a mathematical bifurcation is a qualitative change of asymptotic dynamics 
due to a smooth and slow change in the values of the parameters. A typical example is the famous 
Hopf bifurcation that describes a change from an equilibrium to oscillations. These qualitative 
changes can be continuous and slow (bifurcation theory) but can be also fast and even very fast 
(catastrophe theory). 

Let us recall that “bifurcation theory studies and classifies phenomena characterized by sudden 
shifts in behavior arising from small changes in circumstances, analyzing how the qualitative nature 
of equation solutions depends on the parameters that appear in the equation” ([|I]|). Analysis of 
complex bifurcation points is a classical subject of the theory of dynamical system tlfTll. 

Mathematical Catastrophe theory was invented by Rene Thom in the book Stabilite Structurelle 
et Morphogenese, published in 1972 lfT4]| . He proposed a classification of simple catastrophes 
using geometric properties of surfaces in multidimensional spaces. By E.C. Zeeman “the world 
is full of sudden transformations and unpredictable divergences” that are subjects in Catastrophe 
theory [|20ll . One of the possible explanations of such “sudden transformations” is the existence of 
fast sub-processes into an original system, that relax very quickly to its new asymptotic state. This 
fast relaxation is very hard to observe. 

A set of parameters where a bifurcation or a catastrophe occurs is called a bifurcation set B. 
It is a subset of the space of parameters P and the dimension of the bifurcation set is less than 
the dimension of P. Typically it is a surface in P. The dimension of P (dimP) is the number of 
independent parameters. Because dim B is less then dim P, the number of independent parameters 
of B is less than dim P, and a corresponding part of the parameters m = dim P — dim B depends 
on others. 

We call this phenomenon as “bifurcational parametric simplification 
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The number m is ealled a bifurcational simplification number. 

Conjecture 1. The maximal bifurcational simplification number corresponds to a bifurcation of 
“maximal complexity.” 

The eonjeeture is more philosophieal than it looks. The maximal bifureational simplifieation 
number is elear, but what is a bifurcation of maximal complexity? This question arises beeause 
there exist different types of bifureations: ehange of stability, Hopf bifureation, multiplieity of 
limit eyeles, non loeal bifureations ete... In our opinion, a more aeeurate understanding of the 
maximal bifurcation complexity depends on modeling phenomena. 

The bifureational parametrie simplifieation was first deseribed in papers of Yablonsky et al. 
ifTSll . ifT^ under the name of “eritieal simplifieation” for the famous Langmuir absorption model 
of heterogeneous eatalytie kineties. They observed an essential simplifieation of the kineties at 
extinetion and ignition points that permitted them to find simple eonneetions between the values of 
“bifureation parameters” and kinetie eonstants IfT^ . lfT9ll . As “bifureation parameters” were used 
the maximal kinetie rate i? 3 , is an observable quantity. At ignition and extinetion points is equal 
to the kinetie eonstants of absorption and desorption, that in turn ean be employed to applieations 
involving kinetie eonstant measurements. 

The Langmuir absorption model eombine fast and slow kinetie proeesses. Ignition and extine- 
tion points are points of fast bifureations (eatastrophes). 

1.1. Discussion about ’’maximal complexity” 

We start from an historieal remark. In 1944, Lev Landau notieed that near the loss of stability the 
amplitude of the emergent ’’prineipal motion” satisfies a very simple equation O. It is an example 
of the ’’bifureational parametrie simplifieation”. 

In the ease of an isothermie detailed kineties a eorresponding bifureation set for steady states 
multiplieity is an algebraie variety. By Whitney theorem ETI any algebraie variety admits a Whit¬ 
ney stratifieation. It is natural to eonjeeture that the stratum of minimal dimension eorresponds 
to the ’’maximal eomplexity”. Perhaps numerieal algorithms for the Whitney stratifieation will be 
useful for an evaluation of the ’’maximal eomplexity”. It seems to us that sueh algorithms have to 
be eombined with a subdivision to fast and slow sub-proeesses. 

In the applied bifureation theory there are eonstrueted some interesting algorithms for numer- 
ieal evaluation of the stratum of minimal dimension (maximal eo-dimension), see ffT^ and [fTOll 
that ean be useful to an evaluation of the ’’maximal eomplexity” and deseription of many possible 
transition trajeetories in the vieinity of eritieal points. It was done for isothermal ehemieal system 
El. Of eourse, sueh analysis of sophistieated near-steady-state behavior eombined with a separa¬ 
tion of fast and slow sub-proeesses looks natural. However it still does not provide with a robust 
knowledge on relationships between model parameters whieh determine this eritieal behavior. 

1.2. Bifurcational parametric simplification of chemical kinetic models 

We propose here an aeeurate mathematieal deseription of the bifureational parametrie simplifiea¬ 
tion for kinetie models in the ease of the simplest possible bifureation of equalibriums (steady 
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state) multiplicity. Our justification is based on the bifurcation theory of dynamic systems. We 
also use slow-fast dynamics, i.e. a coordinate free version of the singular perturbation theory, as a 
technical tool for an accurate evaluation of asymptotic analytic expressions that combine different 
reaction parameters |I51 . 

The following questions play a key role in this study: 

• How do the characteristics of multiple steady states and its bifurcations depend on the kinetic 
parameters? 

• How can process properties be readily related to experimental results? 

Let us explain more formally the main ideas of the bifurcational parametric simplification analysis 
in the case of homogeneous kinetic systems. 

In the vector notation the system of governing equations of a homogeneous kinetic system can 
be written as 

where ^ = (^i,-02, •••, is the thermochemical state vector and K = {ki, k 2 ,kg) is the 
kinetic parameter vector. (A detailed description of such models can be found in the next section.) 

Suppose our model has multiple steady states ■01,s, '02,s, •••, 0z,s that are all consistent solutions 
of the functional system 

O = F(0,iT). (1.1) 

These steady states 02,s, •••, 0«,s are functions of the kinetic parameter vector iT = (/ci, /c 2 ,..., kg), 

i.e. 

0i,s = 02,s = 92 {K), ..., i/ji^g = gi{K). 

The simplest possible bifurcation is a coincidence of two steady states 0i_s = gi{K) = g 2 {K) = 

02,s. These bifurcations can be termed as multiplicity bifurcations. The multiplicity bifurcation 
happens if all coordinates of both steady states coincide, i.e. we obtain n functional/algebraic 
equations 

gi{kx, /c2, kg) = g2{ki, k2 ,..., kg) (1.2) 

between kinetic parameters ki, k 2 , ■■■, kg. Therefore there exists at least one functional dependence 
between the kinetic parameters and at most n functional dependencies. The number of independent 
kinetic constants / depends on the structure of the steady states system [LTl 

For this case Q = n — I is the bifurcational simplification number. Any additional bifurcation 
point produces an additional functional equation for kinetic parameters. 

The bifurcational parametric simplification with the maximal simplification number Q corre¬ 
sponds to a bifurcation of “maximal complexity.” It is clear that in our case the maximum criti¬ 
cal simplification is just a set of kinetic parameters. It is just a set of kinetic constants (a value 
Kq = (fci,o, ^ 2 , 0 ) •••) kgfl)) for which there exists a maximal number of functional dependencies. 

(fL2l) . 
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Realization of the bifurcational parametric simplification analysis of kinetic 

systems 

Here we propose an efficient and quite simple algorithm adapted to models of chemical kinetics. 
Main ingredients of this algorithm are: 

- coordinate transformation that reverts the original model to the slow-fast onellSl: 

- the method of invariant manifolds for slow-fast dynamicsl|4l; 

- bifurcation analysis of invariant manifolds PTl 

We apply this algorithm to the Langmuir mechanism. We corrected some inaccuracies of pre¬ 
vious studies by Yablonsky et al using a more delicate algorithm for revealing slow-fast dynamics 
a, in particular we found that the point (0, 0) is not a steady-state. The main invariant of slow 
invariant manifolds is described as a bifurcation of “maximal complexity” of the Langmuir model. 

Possible practical applications 

A practical application is usually the weak point of any theoretical algorithm. We hope that ex¬ 
perimental studies of practical systems for bifurcational kinetic parameters can help in evaluating 
kinetic parameters by using analytic expressions for observable quantities. These analytical ex¬ 
pressions exist due to the bifurcational parametric simplifications phenomenon . 

General Conclusions and discussion 

• Any simple bifurcation point produces at least one functional dependence between kinetic 
parameters. 

• A complex bifurcation point (for example, coincidence of more then two steady states) pro¬ 
duces a number of functional dependencies that correspond in complexity. 

• At a bifurcation point of “maximal complexity” the number of independent parameters is 
minimal. This means that all other kinetic parameters can be analytically recalculated. 

• It is an open question as to how one can use bifurcation critical simplification for experimen¬ 
tal study, but at least for ignition an extinction points it looks possible. 

Conclusions about tbe Langmuir model 

• All steady states, their type, stability and attraction domains are classified. Results about 
type and stability of steady states mainly confirm previous study ifTSll but are more accurate. 
It was proved that (0, 0) is not a steady state as was claimed in [[T5ll . Attraction domains have 
not been previously studied. 

• All bifurcation points which are points of the multiplicity change have been studied. Simi¬ 
larly to IfTSl . ignition and extinction bifurcation points have been observed. Types of these 
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points were analyzed using the speeial eoordinate transformation that is a subjeet of the 
theory of singular perturbed veetor fields [|5]|. 

• Bifureation parameters B 2 for ignition and Bi for extinetion were introdueed and analyzed. 
At bifureation points, eorresponding parameters are vanishing. Bifureation points are ana¬ 
lyzed in details. Results are illustrated with the help of invariant slow manifolds (eurves) by 
a number of figures. 

• At the bifureation point of “maximal eomplexity” both parameters are vanishing. This point 
provides with maximal eritieal simplifieation. This type of bifureation, whieh previously 
was not observed or analyzed, plays a erueial role in our study. 

• The first approximation of slow invariant manifold is analyzed mainly in a vieinity of the 
point (0,0). This analysis permits us to qualify the bifureation type of “maximal eomplexity” 
and to evaluate reaetion rates. 


2. Langmuir mechanism 

It is well known that the simplest meehanism for interpreting isothermal eritieal phenomena in het¬ 
erogeneous eatalysis is the Langmuir meehanism. It is a typieal non-linear three -step adsorption 
meehanism. For CO oxidation on platinum, this meehanism is written as follows: 

1. 2 Z + O 2 ^ 2ZO; 

2. Z + CO^ ZCO\ 

3. ZO + ZCO ^2Z + CO 2 . 

Here Z is the free aetive eatalyst, ZO and ZCO are speeies adsorbed on the eatalyst surfaee. The 
first and seeond step are eonsidered to be reversible. 

We will divide the analysis of the problem into two parts: linear analysis (analysis of the 
stoiehiometry matrix) and non-linear analysis (asymptotie analysis of bifureations). The dimension 
of a kernel of the stoiehiometrie matrix give the relationship between rates of reaetion; from the 
bifureation analysis we obtain an additional relationship between reaetion rates. 

We study a more simple ease of irreversible first reaetion. It is valid for wide domain of 
eoneentrations and pressures [[T^ 

1. 2 Z + O 2 ^ 2ZO; 

2. Z + CO^ ZCO; 

3. ZO + ZCO ^2Z + CO 2 . 
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2.1. Kernel of the stochiometric matrix 


In this case the stoiehiometrie matrix for x := [ZO]^y := [ZCO]^z := [Z] is 

/ 2 0 0 -1 \ 

^= 0 1 - 1-1 

\ -2 -1 1 2 / 


Combining the stoiehiometry matrix S with the rate reaetion veetor (i?i, R- 2 -i R 3 ) gives us the 

following kinetie model in terms of kinetie rates 


dx 

dt 


— 2-Ri — R 3 


dt 


i?2 — R-2 — R 3 


—— — —2Ri — R 2 + R-2 + 2i?3 
dt 

Here Ri = kiPo^z^, R 2 = k 2 Pcoz, R -2 = k- 2 y, R 3 = k^xy. 

By standard ealeulations the rank of the stoiehiometry matrix S equals two. It means that the 
third row of the matrix is a linear eombination on the first and seeond rows, R 0 W 3 = —{Rowi + 
R 0 W 2 ). This relation gives us the existenee of a well known linear integral whieh relates to eatalyst 
mass eonservation low ifT^ 


X + y + z = 1. 

Beeause rank S is 2, the system has no other linear integrals. 

Henee, the eorresponding kinetie model of the open eatalytie system ean be redueed to 

doc 

— = 2 kiPo^{l -x-yf - k^xy 


dt 


k 2 Pcoi)- -x-y)- k- 2 y - k^xy. 


This system has an invariant domain U = [0,1] x [0,1] in the x, r/-phase plane. 
The eorresponding steady state model is: 


2Ri — i?3 


( 2 . 1 ) 


R 2 — R-2 — -Rs- (2.2) 

The partial pressures P 02 Pco are eonsidered to be parameters of our model. 

Conclusion. Equations (12.11) and (12.21) show that only two of four system parameters kiPo 2 , ^ 2 ^ 00 , k- 2 ,k 3 
are independent at any steady state. 
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We rewrite the steady state system (I2.1[|2.2I) in (x, y) phase coordinates 

2kiPo^{l -x-yf - ksxy = 0 

fc 2 -Pco(l - X -y) - k- 2 y - ksxy = 0 . 

This system has an obvious steady state solution (1, 0). By simple calculations we obtain the 
following cubic equation for other possible steady states: 

-2kiPo2klx^ + {2kiPo^kl - AkiPo2k-2k3 - klk2Pco)x^ 


+ {AkiPo2k-2k3 - 2/ciPo2^-2 “ k^klP^o “ k_2k3k2Pco)x + 2fciPo2^-2 = 0- 

This is a cubic equation for the variable x. Corresponding values of y can be evaluated from the 
previous system. 

In typical catalytic processes the elementary reaction between two species on the catalytic 
surface x and y {ZO + ZCO 2Z + CO 2 ) is fast, i.e. ks » /ciPoa, k 2 Pco, k- 2 . This permits 
us to introduce a small parameter e := l/k^ and write the discriminant of the previous cubic 
equation in the following form: 

(—2/C1P02 + k2PcoY{k‘^-2 ~ 8kiPo2k-2 + 2k2Pcok-2 + k^P^o) — 


eAkiPo2{2Ak\PQ^k^_2 + k2Pco{k2Pco + k-2){2k\P^o + 2/c2-Pco^-2 ~ ^- 2 )” 
2/ClPo2^-2(10fc2Pc(9 + k2Pcok-2 + k\)) 


—s‘^AkiPQ^k2Pcok-2{~k‘^2 + “ 1 “ 8k-2{8kiPo2 + k2Pco)) 


t co'^- 2 - 

In the zero approximation £ = 0, the discriminant becomes very simple 


{k2Pco ~ ‘2kiPo2Y{k‘^2 ~ 8kiPo2k-2 + ‘2k2Pcok-2 + k^Pco)- ( 2 - 3 ) 


It can be presented as a product of two terms: ^ 2 Pco—2ki P 02 and k '^2 ~ 8^1 -P 02 A:_ 2 +2^ 2 Peak -2 + 
k^Pco- The corresponding quantities 


P2 := 


Pi := 1 - 


k2Pco . 

2kiPo2 ’ 


fc ^2 - 8kiPo2k-2 + 2k2Pcok-2 + = 

= {k2Pco + k-2f - 8kiPo2k-2, 


(2.4) 


will be used later as bifurcation parameters. A bifurcation happens if Pi = 0, or P 2 = 0. If 
Pi = P 2 = 0 , then the corresponding bifurcation has maximal complexity. 
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3. Transformation of the Langmuir system to a slow-fast sys¬ 
tem 

The fast reaction rate k^xy is presented in both equations. The original system does not have any 
fast-slow division in the original coordinate system {x, y). 

By the main concept in the theory of singularly perturbed vector fields (Sll, BUl, it is possible 
to construct another orthogonal coordinate system where the original model becomes a slow-fast 
one. In the case of the Langmuir system it is only a simple rotation. 

An obvious candidate to a small parameter is e := 1 /^ 3 . We use the following orthogonal 
coordinates transformation 

X — y 

X + y 
V = — 

V2 

With respect to the new coordinates, the system has the form 

^ = \/2/ciPo2(1 - - ^v) + - u), 

^ ^ (^AkiPo^il - \/2uf + 2k2PcoO- ~ - V2k-2{v - u)^ . 

This is a standard slow-fast system (singularly perturbed system) that combines singular and reg¬ 
ular perturbations in the second equation. The singular perturbation induces the small parameter £ 
before and the regular perturbation induces the same small parameter before the second term 
in the right hand side of the second equation. 

3.1. Zero approximation £ = 0 for the singular perturbation 

At this stage, all the machinery in the standard singular perturbation theory can be applied to the 
analysis of the system, including its decomposition to slow and fast subsystems and reduction of 
the near steady state dynamics to dynamics on invariant slow manifolds (slow curves in our model). 

The zero approximation e = 0 (k^ ^ 00 ) of slow invariant curve is given by equating e to zero 
before ^ at the second equation. The slow curves equation is 


To simplify analysis of the steady states we put e = 0 for the regular perturbation in this 
equation. The corresponding approximation of slow curves is represented by = 0. It has 

two branches u + v = 0 and u — v = 0, or in the original coordinates a; = 0 and y = 0. Both 
branches of the slow invariant manifold are stable (attractive). 


2^/2 


{AkiPo,{l - V2uf + 2 k 2 Pco{l - V2v) - V 2 k_ 2 {v - «)) = 0 . (3.1) 
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Let us remark that the two branehes u + u = 0 and u — v = 0 have an intersection point (0, 0) 
where the last approximation is problematic. Around the point (0, 0) a more accurate approxima¬ 
tion is necessary. It is easy to check that (0,0) is not a steady state of the original system (se d2.1[ 
12.21) . This problem is discussed in detail in Appendix 2. 

In the original coordinates the branches are a: = 0 and y = 0. Hence any steady state at this 
approximation has zero of its coordinates, and can not be used for the evaluation of the reaction rate 
= k^xy at any steady state, because it gives the reaction rate zero. This means that we have to 
use a more accurate approximation 13 .1 1 for = k^xy. We shall use the steady state relationships 

of reaction rates (12.1112.21) for evaluations of R 3 . 


3.2. Multiplicity of bifurcation parameters 

I) For the first branch v = uor y = 0 of the slow curve the multiplicity bifurcation parameter is 




k2Pco 

‘^kiPo2 


('see 12.41) . 

There are only one or two steady states on the branch v 
system. The first one is (1,0) and does not depends on Bi. 
The second one is 



k2Pco 

2kiPo2 



uor y 


0 in the original coordinate 


If i?i > 0 this steady state is stable. If i?i < 0 then the second steady state does not belong to the 
invariant domain U, i.e it has no physical meaning. The bifurcation value of Bi is equal 0, i.e in 
the bifurcation point 

k 2 Pco = ‘^kiPo 2 - (3-2) 


Conclusion At the bifurcation point Bi = 0, the parameters kiPo 2 k 2 Pco dependent. 
This dependence is a bifurcational simplification. The change of sign Bi from positive to negative 
corresponds to extinction, that follows from the evaluation of the reaction rates (next subsection). 

II) For the second branch v = —u or a: = 0 of the slow curve the multiplicity bifurcation 
parameter is 

B 2 = (/c- 2 )^ + 2fc_2 {k2Pco ~ ^kiPo2) + {k2PcoY 

(see l2.4l) . 

Let us remark that parameters Bi and B 2 are not completely independent. By simple calcula¬ 
tion we obtain that B 2 = (/c- 2 )^ + 2 /c _2 (- 4 /ciPo 2 -Bi - k 2 Pco) + {k 2 Pcof = {k 2 Pco - k- 2 f - 
8 /c_ 2 /ciPo 2 -Si- From this relation we conclude that if Bi < 0,then P 2 > 0 and if i?i = 0 we have 
B 2 = {k 2 Pco — k- 2 )‘^ > 0 , and the situation Pi < 0 , P 2 < 0 is not possible. 

There exist two, one or zero steady states on this branch, depending on the sign of B 2 . For 
P 2 < 0 steady states does not exist. Existence of one steady states corresponds to P 2 = 0. For 
P 2 > 0 there exist two steady states on a; = 0. Let us remark that /c _2 is the leading ki netic 
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parameter, beeause for /c _2 = 0 the bifureation parameter B 2 is always positive and the bifureation 
does not exists. 

Conclusion. At the bifureation point i ?2 = 0 there exists a funetional relation between param¬ 
eters kiPo 2 , fc- 2 - It is a bifureational simplifieation as well. The ehange in the sign of B 2 

from positive to negative eorresponds to ignition, that follows from the evaluation of the reaetion 
rates (next subseetion). 


3.3. Calculation of reaction rates 

We ealeulate the reaetion rates at the steady states on the first slow eurve y = 0 using the zero 
approximation for Ri and R 2 and the steady state relations between rates R 3 = 2 Ri and Rs = 
R2 — R-2 for the evaluations of R^ and R_2 . 

We use the following standard representation of the steady states coordinate Xg, Ps based on 
the regular perturbation theory: 

{xs = xq + exx -f ..., 2 /s = 2/0 + ^2/1 + •••)• 

Hence R^ = ^((xq + exi + ...)( 2 /o + ^ 2 /i + •••))• approximation e = Qwe have xq = 0 or 

2/0 = 0. Therefore R3 = xipo or R3 = xopi. 


The first branch of the slow curve is 2 / = 0. 

At the steady state (1,0) all reaetion rates equal zero. The seeond steady state 


k 2 Pco 


,0 


2kiPo2 

exists for Bi > 0. At this steady state the reaetion rates are 


Ri = 


R 2 = 


klP^o 
AkiPo 2 

klPho 

2 kiPo 2 

R-2 = 0 , 


Clearly, we have an equality 


Rz 


2 Ri = 


U2 p2 

2kiPo2 


2Ri — R 2 — R-2 


In this approximation the rate of the seeond reverse reaetion at the first braneh is zero. Henee all 
reaetions of the detailed meehanism ean be eonsidered to be irreversible. 
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Remark. In the first regular approximation (xq + exi), at the first braneh the rate of CO des¬ 
orption (the seeond reverse reaction) is zero. But for more accurate approximations it is positive, 
but small enough. Our approximation is of the order 0(e) for singular perturbations. 

The second branch of the slow curve is x = 0. 

At the case i ?2 > 0. There exist two steady states (on this branch): 

/c 2 + k2PcO i 


0 , 1 - 


4fciP, 


02 


Corresponding reaction rates are: 


1 TD [ ^-2 + k2PcO i s/B2\ 

- ) ■ 


Ry — koP( 


2 -rcO 


fc 2 + k2PcO ± 




R—2 — k—2 


= 2Pi = 2fciP, 


AkiPo2 J 

k_2 + k2Pco 4 = \/^\ 


O 2 


4fciPo2 

k 2 + k2PcO i 






4fciPo2 ) 

The equation 2Pi = R 2 — R -2 is also valid for the second branch of the slow curve. 

Bifurcation points. 

Now we will calculate the steady-state reaction rates at the bifurcation points. 

At the first bifurcation point Pi = 0 (extinction point) the reaction rates for the branch y = 0 


are: 


Pi = /C 1 P 02 , 


P2 — 2/C1P02, 


P-2 — 0, 


P 3 = 2Pi = 2kiPo2- 

After this bifurcation the system ’’jumps” to the stable steady state on the branch x = 0. The 
condition Pi = 0 is fulfilled at the reaction rate 2 kiPo 2 j the extinction 

point. Thus, the main reaction rate of P 3 at this steady state is 


-^3,1 — 


{k-2 


2kiPc 


02 
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Because ’’jumps” from the larger value 2 kiPo 2 to the smaller value 2 kiPo bifurcation 
point is the extinction point, i.e it ’’jumps” from the extinction value Ra to the ’’after extinction” 
value By previous calculations 


RaRb — (^- 2 )^ 


(3.3) 


At the second bifurcation point i ?2 = 0 (i.e. {k 2 Pco + ^- 2 )^ = 8 k_ 2 kiPo 2 ) reaction rates for 
the branch x = 0 are: 

/C_2 

Ri = —, 

2 


R2 = 


k2PcO n , 7 T-, \ k2PcO 


AkiPo 2 


ik.2 + k2Pco2) = 


4A:iP( 


02 


R -2 = {Ak,Po 2 - {k -2 + k 2 Pco)) = 


^hPc 


02 


^kiP, 


02 


AkiPo2 ~ V8kiPo2k-‘^ , 


R 3 — 2i?i — k- 2 - 

After this bifurcation the system ’’jumps” to the stable steady state on the branch y = 0. The 
main reaction rate P 3 at this steady state is 

k^P^Q 

— 


2 kiPc 


02 


^ 2 Eco_ this bifurcation 


Because P 3 ’’jumps” from the smaller value k -2 to the bigger value „ 

ZfZ\ XQ2 

point is the ignition point, i.e it ’’jumps” from the ignition value Pc to the after ignition value Rd- 
By previous calculations 


k- 2 klP^o 

At the bifurcation point of the maximal complexity Bi = B 2 = Owe have a further simplifica¬ 
tion: 

R 3 = k-2 = k2Pco = 2/ciPo2- 


4. Dynamics of Langmuir system 

In all figures, straight lines represent the zero approximation e = 0, dashed lines represent the first 
approximation, thick points are steady states, and small arrows represent system trajectories. All 
figures are produced by the computer program “Mathematica”. 
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By standard stability analysis we have that for Bi > 0 and B 2 > 0 the steady states of (1,0) 

and 

/p, . k-2 + ^2-fc'O + 

V ’ J 

are saddle points. The steady states of oj and 



k_2 + k2Pco ~ y/^\ 
AkiPo^ ) 


are stable (attractive) nodes. 

If the initial conditions satisfy the following inequality: 


X < y — 1 + 


k-2 + k2Pco + \/ B 2 

AkiPo2 


all trajectories are attracted to the node 



k-2 + k2Pco — \/^\ 
4 /C 1 P 02 / 


For any other initial conditions all trajectories are attracted to the node 


(see Figure 1). 


k 2 Pco 
2fciPo2’ J 



Figure 1: Dynamics of the system in the case Pi > 0, P 2 > 0. Solid line: zero approximation of 
slow manifold, dotted line: first approximation of slow manifold, red points: steady states. 


In the case Pi > 0, P 2 < 0 there are two steady states on the branch y = 0 and no steady state 
on the branch x = 0. Hence, all trajectories are attracted to the node 


1 - 


k 2 Pco 

2 A:iPo,’ J- 
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0.8 

0 . 6 

0.4 

0.2 

0 


Figure 2: Dynamics of the system in the case Bi > 0, i ?2 < 0. Solid line: zero approximation of 
slow manifold, dotted line: first approximation of slow manifold, red points: steady states. 

(see Figure 2). 

In the case i?i < 0, i ?2 > 0 the steady state is 

k-2 + k2PcO + V -^2^ 

4fciPo3 ) ’ 

and it does not belong to the invariant domain U. Hence we have only two singular points. 

The first one (1, 0) is a saddle point, while the second one is a stable node 

/c_2 + /c2-fcO ~ \fB2 
AkiPo2 

All trajectories are attracted to the second point (see Figure 3). 



Figure 3: Dynamics of the system in the case i?i < 0, i ?2 > 0. Solid line: zero approximation of 
slow manifold, dotted line: first approximation of slow manifold, red points: steady states. 

In bifurcation cases where one or both parameters Bi, B 2 become zero more accurate approx¬ 
imations are necessary and such analyses are beyond the scope of the zero approximation. 
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Conclusion. If i ?2 > 0 and the sign of the bifurcation parameter Bi is changed from the 
positive to negative one, the reaction rate decreases dramatically from large to small values 
(compare figures 1 and 3). It means that Bi = 0 corresponds to the extinction point. 

If i?i > 0 and the sign of the bifurcation parameter B 2 is changed from the positive to negative 
one, the reaction rate R 3 make a ’’jump” from small to larg values (compare figures 2 and 3). This 
means that B 2 = 0 corresponds to the ignition point. 

It is shown that the point of maximal bifurcation complexity present a coincidence of ignition 
and extinction points.” 


5. About the phisyco-chemical meaning of bifurcation param¬ 
eters and bifurcations conditions 


We follow chapter 9 of Marin and Yablonsky’s book lfT9ll . 

The condition i?i = 0 is equivalent to z = 2 kiPo ~ ^ corresponding equilibrium, 

i.e. the entire catalyst surface is empty at the equilibrium. Therefore, the bifurcation parameter 
Bi ■= 1 — 2 kiPo determines the difference between the complete concentration of catalytic sites, 
i.e,unity, and the corresponding steady-state concentration of the empty catalytic sites. 

The condition i ?2 = 0 can be reformulated in terms of the parameter 


{hPco + fc-2) 
8kiPo2 


The condition i ?2 = 0 is equivalent to H = 1. 

By the previous analysis the condition Bi = 0 determines the extinction point and the condition 
^2 = 0 or TT = 1 determines the ignition point. 

Detailed explanations of the condition H = 1 can be found in the book ( |[T9ll , chapter 9). We 
repeat here these arguments. “At the ignition point, the steady state reaction rate R^ = k^xy is 
determined only by the CO desorption rate coefficient. It does not depend on the composition of 
the gas-phase mixture.” At the ignition point the concentration Zig of active sites is equivalent to 
2zeg, where Zeq is an equilibrium concentration of free active sites. Let us discuss in more detail 
the concept of equilibrium concentration. Suppose that in our system we have only one reversible 
adsorption step, i.e. an interaction between the empty catalytic sites and gaseous CO. Then at the 
equilibrium which in this case is the steady state as well 


k2Pco^eq 2l/eg- 


From these equations, obviously 

_ ^-2 _ 1 

“ k2PcO + k-2 ~ Keg,2PcO + l’ 
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where A"eq ,2 is the equilibrium eonstant of CO adsorption. 

Before ignition, the steady state concentration of adsorbed oxygen, x is very low, therefore in 
the vicinity of ignition we approximately have 

y 1 ^ign 1 ^Z. 

At the point of the maximal bifurcational complexity both conditions i?i = 0 and if = 1 are 
fulfilled. It means that at this point two phenomena, ignition and extinction, merge. All the 
catalyst sites are empty. Using z = 2zeq = 1, we have 

_ 2 _ 

“ K^.2Pco + 1 “ ^ 


and Keq^ 2 Pco = 1- 

The same equation can be obtained in a more similar way, because at the point of the maxi¬ 
mal bifurcational complexity (MBC-point) ignition and extinction are merged, and the following 
equation is valid 

Rext Rign 

where Rext and Rign are reaction rates at ignition and extinction points respectively. Consequently 

k -2 = hPco and Keq, 2 Pco = 1- 

The last equation is quite interesting. It is nothing but the ’degenerate’ Langmuir equation in 
the limit case when all the catalytic surface is covered by CO (rather to say almost covered). 

This equation gives a unique possibility to estimate Keq based just on the MBC -point exper¬ 
imentally observed: 

Reg,2 = PcO,MBC, 

where Pco,mbc is the partial pressure of CO at the MBC-point 

Remark Distinguishing the well-observed and ill observed critical characteristics and param¬ 
eters. 

In catalytic CO oxidation over the Pt catalyst, the reaction rates at ignition and extinction 
points are ill -observed, but the after-ignition and after-extinction points are well-observed. It is 
possible to find reaction rates of the ill-defined points (ignition and extinction) and parameters 
of adsorption steps based on the information regarding the well-defined points (after-ignition and 
after-extinction). 


6. Conclusions 

A phenomenon of the bifurcational parametric simplification was presented as a generalization 
and a justification of the critical simplification principle described in [fT^ . The parametric sim¬ 
plification was applied to the Langmuir model using its transformation to the slow-fast dynamic 
system based on the concept of singularly perturbed vector fields |l5l. A detailed analysis of the 
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Langmuir model was performed, and a bifureation of the maximum eomplexity was presented as 
a new peeuliarity in terms of kinetic parameters. 


Appendix 1. Mathematical models and the structure of chemical 
reaction mechanisms 

Additional definitions of the chemical source term are required for the exposition. A pure homoge¬ 
neous system is considered. It is represented by a system of ordinary differential equations (ODEs) 
that describes the mechanism of chemical kinetics by a system based on the mass action law. 

The mathematical model describes the temporal evolution of chemical state vector, where ^|Jj 
represents a concentration of the its chemical substance. A dimension of the system is the number 
of species (reactants), n: 

In vector notation the system of governing equations of a homogeneous system can be written 
in autonomous form as (see e.g. mEm) 

^ e C (6.1) 

Here the so called chemical source F represents the chemical mechanism: reactants partic¬ 
ipate in Ur elementary chemical reactions. Pressure and enthalpy are considered to be constant. 
The elementary reactions are 

i = 1, ...,nr : ai^iAi -\- ... -\- ans,iAn^ -H- + ... -f (3ns,iAn^ (6-2) 

where Ai,,, chemical species which participate in Ur elementary reactions. Reaction rates of 
these reactions relate to the mass action law which was already mentioned [7]. This law implies 
the polynomial form 


Tis ris 

Ri{ip) := kf JJ - K JJ i = 1, ...,nr. (6.3) 

i=i i=i 

The rate constants kf, k~ are characterized by exponential Arrhenius dependences on the temper¬ 
ature. 

Therefore, components of F{ip) are composed in the following way 


Tlr 

= (^)) (6-4) 

i=l 

where := (3j^i — aj^i are components of the so-called stoichiometric coefficient (n* x n^) matrix 
S, and i?(^) is a non-linear vector function of the elementary reaction rates. 
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A specific structure of models of chemical kinetics is discussed now. The structure is based on 
the stoichiometric matrix S. It is typical in chemical kinetics that the number of reactions is much 
larger then the number of species i.e. rir S> rig and S' is a rectangular matrix of a linear mapping 
from the reaction space to the system state space . The vector of elementary reaction rates 
is a non-linear map from the state space i?”" to the reaction space i?”’’. Their composition 
represents the chemical source term (vector) F(\E^) that maps the state space i?”® into itself. 

This special structure of kinetic models (16.11) - (16.41) can be illustrated by the following diagram 


RUr 

m) 

^ I s 

R"’ —> R"' 

Fm 


The initial linear reduction procedure is possible for the linear map represented by the matrix 
S. Denote by 


kei S = {w e R^^\S{w) = 0} 


the kernel of the linear map S : i?”’’ —)■ i?”®. It is clear that dim(ker S) > — n^. If dim(ker S) > 

Ur — rig, then there exist 

q = dim(ker S) — {jir — rig) > 0 


linear integrals in the original system. Hence, without loss of generality, the first rig — q lines of 
S := ||7yi|| can be assumed to be linearly independent, while the q lines k = rig — q + 1^ rig are 
linearly dependent 

Us-q 

Ik,- ^ ^ 

OL = l 

As the result, the following linear integrals of the original system can be found, given by the 
(g X rig) matrix 



/ cri,l 

g 1 0 

0 \ 

Zc = 

<^2,1 

^2,ns—q 0 1 

0 


\ ^9,1 

^q,ns—q 0 0 

-1 / 



o' 

III 

o 


which means Zc defines the conserved quantities of the system 


Zc 

d^! 

dt 

Z, • = Z,-S- R{ir) 

= 0, 

i.e. 


1 

O 

O 
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This procedure permits us to reduee the dimension of the system using the linear integrals whieh 
typieally have a simple physie-ehemieal meaning eorresponding to eonserved quantities of ehemi- 
eal elements. 

A detailed study of ker S has a long history. It is known as the “methabolie pathway analysis 
of null spaee of stoehiometrie matrix” (see, for example, lIHl). A more detailed analysis of linear- 
algebraie approaehes with help of moleeular and Horiuti matriees ean be found in IfTTlI . 

Steady states ...^m of this model are solutions of the following analytie equation 

F('l/i) = 0,i = 1,2, ...,m. 

Its eoordinates are funetions of main system parameters. Therefore, the linear integrals produee 
funetional dependeneies between main model parameters. Solving these funetional dependeneies 
we reduee number of main independent parameters by the number of the linear integrals. Typieally 
explieit solutions of the steady state equations and its eonsequenees for main system parameters, 
have a very eomplieated analytie nature beeause of the high non-linearity of the original models. 

It is a first step of the parametrie eritieal simplifieation. 

The seeond step of the parametrie eritieal simplifieation is eonneeted with bifureations. The 
notion of a bifureation is not eompletely formal. In the ease of steady state multiplieity, it means 
eoineidenee of two different steady states i.e. = ^j,i ^ j- It gives an additional funetional 
dependenee between the main parameters of the model and therefore permits us to reduee at least 
one independent parameter. Beeause bifureation typieally has a simple physieal meaning (igni¬ 
tion or extinetion points, appearanee of oseillations ete.) it permits us to eheek values of leading 
parameters experimentally. 


Appendix 2. About the first approximation of the slow invariant 
curve 

Reeall that in u, v-eoordinates the Langmuir system has the form 

^ = V2hPo,{l - V2v) - ^Pco{l - V2vf, 

^ ^ ^ ^ (4fciPo2(l - + 2k2Pco{^ - . 

Beeause the zero approximation e = 0 is not informative at the interseetion point (0, 0) of the 
two branehes u = v and u = —v, we ealeulate the first approximation of the slow eurve. It is 

-v'^ + e (k 2 Pco{^ - ^ -f 2 /ciPo 2 (l “ = 0 

\ V J 

(see Figure 4). 
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We give here a short sketch of the first approximation of steady state multiplicity. We do 
not present here all of the long elementary calculations, but focus rather on the vicinity of the 
problematic intersection point (0, 0). Far from (0, 0) the first approximation loses its corrective 
influence on the results. 

The first bifurcation parameter Bi in the first approximation becomes 


^ 1,1 = (1 


^2-PcO x2 _ ‘^^2^00 
2kiPo2 kiPo2 


Depending on this bifurcation parameter the system has four, three or two steady states. Existence 
of three steady states ((0,1), (1,0) and 


,1 

'2 



9^2 p2 

■^^ 2 -' CO 

kiPo 2 


1 / jfcjpgo 
2\l kiPo2 


corresponds to = 0. The third steady state 



9^2 p2 

^^ 2 ^CO 


kiP, 


02 


1 j 2klP^co 

2 \l k,Po 2 


belongs to an e-vicinity of (0, 0). For e —)■ 0 this point is disappears from the vicinity of the (0,0) 
point. This steady state disappears for Pi i < 0. 

The first order approximation of the slow manifold is presented at Figure 4. 


V 



Figure 4: Solid line: zero approximation of slow manifold, dotted line: first approximation of slow 
manifold 

In the invariant domain we have that /c 2 Pco(l ~ + 2 kiPo2Y ~ > 0. This 

inequality leads to the conclusion that non trivial steady states lie inside the invariant domain. In 
the case Pi i =0 two non trivial steady states come together and we obtain one steady state that it 
is a saddle-node. In this case, if the initial condition satisfies the inequality: x < y, all trajectories 
attract to the (0,1) point. If the initial condition satisfies the inequality x > y, all trajectories 

attract to the Point. 
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The case of i?i = 0, 52 = 0 is more complex and the standard asymptotic analysis is not 
applicable. Actually, this case is more interesting because this bifurcation point is the point of the 
maximal simplification of the system. At this point we obtain that k _2 = k 2 Pco = 2fci5o2- 
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